-------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
 opened on:  21 Jan 2023, 17:18:51

. log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
 paused on:  21 Jan 2023, 17:18:51
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
resumed on:  21 Jan 2023, 17:18:51

.       sum `demos' ;

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
        male |      1,896    .4641303    .4952851          0          1
         age |      1,896    51.10127    16.33085         18         90
     college |      1,896    .7480163    .4313104          0          1
       white |      1,896    .7675108    .4180058          0          1

.    log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
 paused on:  21 Jan 2023, 17:18:51
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
resumed on:  21 Jan 2023, 17:18:51

. foreach y of varlist `demos' { ;
  2.    replace `y' = 0 if dummy == 0 ;
  3. } ;
(3,576 real changes made)
(7,584 real changes made)
(5,700 real changes made)
(5,860 real changes made)

. log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
 paused on:  21 Jan 2023, 17:18:51
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
-------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitivities.log
  log type:  text
resumed on:  21 Jan 2023, 17:18:51

. eststo clear ;

. /* PREFERENCE SPACE FIXED COEFFICIENT MODEL */
> 
> clogit ce cost ec_g ec_p hr_g hr_p md_g md_p dummy `demos' [pw=wgt_final], group(coid) cluster(resp_id) ;

Iteration 0:   log pseudolikelihood = -3745.4716  
Iteration 1:   log pseudolikelihood =  -3713.491  
Iteration 2:   log pseudolikelihood = -3713.4229  
Iteration 3:   log pseudolikelihood = -3713.4229  

Conditional (fixed-effects) logistic regression

                                                Number of obs     =     15,168
                                                Wald chi2(12)     =     573.77
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -3713.4229               Pseudo R2         =     0.1281

                            (Std. Err. adjusted for 1,896 clusters in resp_id)
------------------------------------------------------------------------------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        cost |  -.0560818   .0032908   -17.04   0.000    -.0625315    -.049632
        ec_g |   .0151355   .0034873     4.34   0.000     .0083005    .0219705
        ec_p |  -.0153156   .0037155    -4.12   0.000    -.0225979   -.0080332
        hr_g |   .0107416   .0028788     3.73   0.000     .0050993    .0163838
        hr_p |  -.0246492   .0084162    -2.93   0.003    -.0411446   -.0081539
        md_g |   .0077462   .0021081     3.67   0.000     .0036144    .0118781
        md_p |  -.0038339   .0045866    -0.84   0.403    -.0128234    .0051557
       dummy |   1.254328   .2183657     5.74   0.000     .8263388    1.682317
        male |  -.1003348   .0981777    -1.02   0.307    -.2927597      .09209
         age |  -.0154611   .0027465    -5.63   0.000    -.0208442    -.010078
     college |   .2852263   .1006662     2.83   0.005     .0879242    .4825284
       white |  -.2640639    .118147    -2.24   0.025    -.4956278   -.0325001
------------------------------------------------------------------------------

. /* #1 - WTP SPACE FIXED COEFFICIENT MODEL */
> 
> mixlogitwtp ce dummy ec_g ec_p hr_g hr_p md_g md_p `demos' [pw=wgt_final], group(coid) id(resp_id) price(mcost) nrep(3) vce(robust) ;

Iteration 0:   log pseudolikelihood =  -3713.234  (not concave)
Iteration 1:   log pseudolikelihood = -3710.9279  (backed up)
Iteration 2:   log pseudolikelihood = -3707.8265  
Iteration 3:   log pseudolikelihood =  -3707.799  
Iteration 4:   log pseudolikelihood =  -3707.799  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    1907.40
Log pseudolikelihood =  -3707.799               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |     24.398   4.108814     5.94   0.000     16.34487  
>   32.45112
        ec_g |     .30217   .0749335     4.03   0.000     .1553031  
>    .449037
        ec_p |  -.3322179   .0898432    -3.70   0.000    -.5083073  
>  -.1561284
        hr_g |    .201323   .0594592     3.39   0.001     .0847851  
>   .3178609
        hr_p |  -.5544499   .1780947    -3.11   0.002    -.9035091  
>  -.2053908
        md_g |   .1289249   .0423051     3.05   0.002     .0460083  
>   .2118414
        md_p |  -.0257647   .0971089    -0.27   0.791    -.2160946  
>   .1645653
        male |  -2.345356   2.206538    -1.06   0.288    -6.670091  
>   1.979378
         age |   -.317045   .0622715    -5.09   0.000     -.439095  
>  -.1949951
     college |    6.69105   2.473587     2.70   0.007     1.842908  
>   11.53919
       white |  -6.462618   3.385864    -1.91   0.056    -13.09879  
>   .1735528
       mcost |  -2.980282    .096514   -30.88   0.000    -3.169445  
>  -2.791118
-------------+------------------------------------------------------
> ----------
SD           |
       mcost |   .5413913   .1654622     3.27   0.001     .2170913  
>   .8656914
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive

. matrix define b0 = e(b) ;

. constraint 1 [SD]_b[mcost] = 0 ;

. matrix list b0 ;

b0[1,13]
          Mean:       Mean:       Mean:       Mean:       Mean:
         dummy        ec_g        ec_p        hr_g        hr_p
y1   24.397995   .30217002  -.33221787   .20132301  -.55444992

          Mean:       Mean:       Mean:       Mean:       Mean:
          md_g        md_p        male         age     college
y1   .12892487  -.02576466  -2.3453562  -.31704504     6.69105

          Mean:       Mean:         SD:
         white       mcost       mcost
y1  -6.4626182  -2.9802815   .54139134

. matrix b0[1,13] = 0 ;

. eststo: mixlogitwtp ce dummy ec_g ec_p hr_g hr_p md_g md_p `demos'
>  [pw=wgt_final], group(coid) id(resp_id) price(mcost) nrep(3) vce(
> robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -3717.9224  
Iteration 1:   log pseudolikelihood = -3713.8661  
Iteration 2:   log pseudolikelihood = -3713.4239  
Iteration 3:   log pseudolikelihood = -3713.4229  
Iteration 4:   log pseudolikelihood = -3713.4229  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    3802.93
Log pseudolikelihood = -3713.4229               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   22.36605   3.458984     6.47   0.000     15.58656  
>   29.14553
        ec_g |   .2698828   .0658171     4.10   0.000     .1408837  
>   .3988819
        ec_p |  -.2730933   .0698794    -3.91   0.000    -.4100545  
>  -.1361322
        hr_g |   .1915341   .0543041     3.53   0.000     .0850999  
>   .2979682
        hr_p |  -.4395233    .154377    -2.85   0.004    -.7420967  
>  -.1369499
        md_g |   .1381239   .0399314     3.46   0.001     .0598598  
>   .2163879
        md_p |   -.068362    .082007    -0.83   0.405    -.2290927  
>   .0923688
        male |  -1.789082   1.747085    -1.02   0.306    -5.213304  
>   1.635141
         age |  -.2756889   .0495696    -5.56   0.000    -.3728435  
>  -.1785343
     college |   5.085902   1.815979     2.80   0.005     1.526649  
>   8.645155
       white |  -4.708553   2.125485    -2.22   0.027    -8.874427  
>  -.5426787
       mcost |  -2.880945   .0586777   -49.10   0.000    -2.995951  
>  -2.765939
-------------+------------------------------------------------------
> ----------
SD           |
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est1 stored)

. local ll1 = e(ll) ;

.  local par1 = e(k) ;

. /* #2 - WTP SPACE RANDOM COEFFICIENT MODEL - BASELINE */
> 
> eststo: mixlogitwtp ce dummy `demos' [pw=wgt_final], group(coid) i
> d(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p md_g md_p) nrep(2
> 50) vce(robust) ;

Iteration 0:   log pseudolikelihood =  -3692.657  (not concave)
Iteration 1:   log pseudolikelihood = -3365.3904  (not concave)
Iteration 2:   log pseudolikelihood = -3359.8195  
Iteration 3:   log pseudolikelihood = -3304.9303  (not concave)
Iteration 4:   log pseudolikelihood = -3290.8903  
Iteration 5:   log pseudolikelihood = -3266.2897  
Iteration 6:   log pseudolikelihood = -3262.7496  
Iteration 7:   log pseudolikelihood = -3262.6003  
Iteration 8:   log pseudolikelihood = -3262.5887  
Iteration 9:   log pseudolikelihood = -3262.5887  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    1087.50
Log pseudolikelihood = -3262.5887               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |    14.9551   3.394468     4.41   0.000     8.302069  
>   21.60814
        male |  -1.630105   1.681037    -0.97   0.332    -4.924877  
>   1.664667
         age |  -.1399781   .0513579    -2.73   0.006    -.2406378  
>  -.0393184
     college |   3.815571   1.700058     2.24   0.025      .483519  
>   7.147622
       white |  -5.957026   1.910814    -3.12   0.002    -9.702152  
>    -2.2119
        ec_g |   .3188457   .0580402     5.49   0.000     .2050889  
>   .4326024
        ec_p |   -.276512   .0586978    -4.71   0.000    -.3915575  
>  -.1614664
        hr_g |   .2616566   .0586321     4.46   0.000     .1467398  
>   .3765734
        hr_p |  -.6814736    .154192    -4.42   0.000    -.9836843  
>   -.379263
        md_g |   .1527001    .037483     4.07   0.000     .0792347  
>   .2261655
        md_p |   -.177933   .0741067    -2.40   0.016    -.3231795  
>  -.0326865
       mcost |  -1.970238    .104847   -18.79   0.000    -2.175734  
>  -1.764742
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   .8297363   .0873132     9.50   0.000     .6586056  
>   1.000867
        ec_p |   .3350668   .1517776     2.21   0.027     .0375881  
>   .6325456
        hr_g |   .5677385   .1131209     5.02   0.000     .3460256  
>   .7894514
        hr_p |  -1.811267   .2630972    -6.88   0.000    -2.326928  
>  -1.295605
        md_g |   .5223477    .103475     5.05   0.000     .3195405  
>   .7251549
        md_p |   1.011984    .109776     9.22   0.000     .7968271  
>   1.227141
       mcost |  -.2035881    .198862    -1.02   0.306    -.5933504  
>   .1861742
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est2 stored)

. local ll2 = e(ll) ;

.  local par2 = e(k) ;

. test ec_g == -ec_p ;

 ( 1)  [Mean]ec_g + [Mean]ec_p = 0

           chi2(  1) =    0.25
         Prob > chi2 =    0.6196

. test hr_g == -hr_p ;

 ( 1)  [Mean]hr_g + [Mean]hr_p = 0

           chi2(  1) =    6.54
         Prob > chi2 =    0.0106

. test md_g == -md_p ;

 ( 1)  [Mean]md_g + [Mean]md_p = 0

           chi2(  1) =    0.09
         Prob > chi2 =    0.7652

. test (ec_g == -ec_p) (hr_g == -hr_p) (md_g == -md_p) ;

 ( 1)  [Mean]ec_g + [Mean]ec_p = 0
 ( 2)  [Mean]hr_g + [Mean]hr_p = 0
 ( 3)  [Mean]md_g + [Mean]md_p = 0

           chi2(  3) =    6.93
         Prob > chi2 =    0.0743

. test (ec_g == -ec_p == hr_g == -hr_p) ;

 ( 1)  [Mean]ec_g + [Mean]ec_p = 0
 ( 2)  [Mean]ec_g - [Mean]hr_g = 0
 ( 3)  [Mean]ec_g + [Mean]hr_p = 0

           chi2(  3) =    6.92
         Prob > chi2 =    0.0745

. test (ec_g == -ec_p == hr_g == -hr_p == md_g == -md_p) ;

 ( 1)  [Mean]ec_g + [Mean]ec_p = 0
 ( 2)  [Mean]ec_g - [Mean]hr_g = 0
 ( 3)  [Mean]ec_g + [Mean]hr_p = 0
 ( 4)  [Mean]ec_g - [Mean]md_g = 0
 ( 5)  [Mean]ec_g + [Mean]md_p = 0

           chi2(  5) =   15.41
         Prob > chi2 =    0.0088

. /* #3 - WTP SPACE RANDOM COEFFICIENT MODEL - NO SAMPLING WEIGHTS I
> N ESTIMATION */
> 
> eststo: mixlogitwtp ce dummy `demos', group(coid) id(resp_id) pric
> e(mcost) rand(ec_g ec_p hr_g hr_p md_g md_p) nrep(250) vce(robust)
>  ;

Iteration 0:   log pseudolikelihood = -4549.3238  (not concave)
Iteration 1:   log pseudolikelihood =  -4154.196  (not concave)
Iteration 2:   log pseudolikelihood = -4115.7077  
Iteration 3:   log pseudolikelihood =  -4054.461  (not concave)
Iteration 4:   log pseudolikelihood = -4042.6229  
Iteration 5:   log pseudolikelihood = -4018.6184  
Iteration 6:   log pseudolikelihood = -4017.5035  
Iteration 7:   log pseudolikelihood = -4017.4358  
Iteration 8:   log pseudolikelihood = -4017.4332  
Iteration 9:   log pseudolikelihood = -4017.4332  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    1078.94
Log pseudolikelihood = -4017.4332               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   12.51853   2.643274     4.74   0.000     7.337811  
>   17.69926
        male |  -1.603326   1.250862    -1.28   0.200    -4.054971  
>   .8483185
         age |   -.088778   .0418333    -2.12   0.034    -.1707698  
>  -.0067861
     college |   3.117144   1.554821     2.00   0.045     .0697495  
>   6.164538
       white |  -4.158524   1.565472    -2.66   0.008    -7.226793  
>  -1.090254
        ec_g |   .3080984   .0402011     7.66   0.000     .2293057  
>    .386891
        ec_p |  -.2231307    .046195    -4.83   0.000    -.3136713  
>  -.1325901
        hr_g |    .278871   .0479645     5.81   0.000     .1848622  
>   .3728797
        hr_p |  -.5335149   .1231046    -4.33   0.000    -.7747954  
>  -.2922344
        md_g |   .1072954     .02759     3.89   0.000     .0532199  
>   .1613709
        md_p |  -.2453533   .0569111    -4.31   0.000    -.3568969  
>  -.1338096
       mcost |  -1.772502   .3051458    -5.81   0.000    -2.370577  
>  -1.174428
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   .7181847   .1254823     5.72   0.000     .4722438  
>   .9641255
        ec_p |   .3178278   .1336433     2.38   0.017     .0558918  
>   .5797637
        hr_g |   .4661308   .1053461     4.42   0.000     .2596562  
>   .6726053
        hr_p |   2.048825    .756086     2.71   0.007      .566924  
>   3.530726
        md_g |   .4764504    .093553     5.09   0.000       .29309  
>   .6598108
        md_p |   .9562349   .0982402     9.73   0.000     .7636877  
>   1.148782
       mcost |   .5160813   .3709358     1.39   0.164    -.2109396  
>   1.243102
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est3 stored)

. local ll3 = e(ll) ;

.  local par3 = e(k) ;

. /* #4 - WTP SPACE RANDOM COEFFICIENT MODEL - ECOSYSTEM CONDITION F
> IRST */
> 
> eststo: mixlogitwtp ce dummy `demos' if block_f == 1 [pw=wgt_final
> ], group(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p m
> d_g md_p) nrep(250) vce(robust) ;

Iteration 0:   log pseudolikelihood = -1854.9715  (not concave)
Iteration 1:   log pseudolikelihood = -1672.2971  (not concave)
Iteration 2:   log pseudolikelihood = -1672.0963  
Iteration 3:   log pseudolikelihood = -1671.6489  (not concave)
Iteration 4:   log pseudolikelihood = -1651.3314  (not concave)
Iteration 5:   log pseudolikelihood = -1640.1896  (not concave)
Iteration 6:   log pseudolikelihood =  -1631.379  (not concave)
Iteration 7:   log pseudolikelihood =  -1628.765  
Iteration 8:   log pseudolikelihood = -1627.2959  
Iteration 9:   log pseudolikelihood = -1627.2498  
Iteration 10:  log pseudolikelihood = -1627.2428  
Iteration 11:  log pseudolikelihood = -1627.2414  
Iteration 12:  log pseudolikelihood = -1627.2414  

Mixed logit model in WTP space                  Number of obs     = 
>      7,672
                                                Wald chi2(12)     = 
>     424.16
Log pseudolikelihood = -1627.2414               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   15.94971   6.744891     2.36   0.018     2.729964  
>   29.16945
        male |  -1.416514   4.179896    -0.34   0.735    -9.608959  
>   6.775931
         age |  -.1281778   .1200142    -1.07   0.286    -.3634013  
>   .1070457
     college |   6.066002    3.87011     1.57   0.117    -1.519273  
>   13.65128
       white |  -8.839291   6.307674    -1.40   0.161     -21.2021  
>   3.523522
        ec_g |   .2935815   .1500554     1.96   0.050    -.0005216  
>   .5876847
        ec_p |  -.2321428   .0899232    -2.58   0.010     -.408389  
>  -.0558966
        hr_g |   .2845286   .0908137     3.13   0.002     .1065369  
>   .4625202
        hr_p |  -.6093764   .3849361    -1.58   0.113    -1.363837  
>   .1450845
        md_g |   .1812847   .0785387     2.31   0.021     .0273518  
>   .3352177
        md_p |  -.0908524    .126686    -0.72   0.473    -.3391524  
>   .1574476
       mcost |  -2.026205   1.042438    -1.94   0.052    -4.069345  
>   .0169353
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   1.099889   .4844809     2.27   0.023     .1503236  
>   2.049454
        ec_p |   .5774249   .2886256     2.00   0.045     .0117291  
>   1.143121
        hr_g |  -.3480591     .46867    -0.74   0.458    -1.266635  
>   .5705173
        hr_p |   1.264014   .7540878     1.68   0.094    -.2139709  
>   2.741999
        md_g |  -.4771369   .2643894    -1.80   0.071    -.9953307  
>   .0410568
        md_p |   1.321971   .8997441     1.47   0.142    -.4414946  
>   3.085438
       mcost |   -.191513   4.192153    -0.05   0.964    -8.407981  
>   8.024955
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est4 stored)

. local ll4 = e(ll) ;

.   local par4 = e(k) ;

. /* #5 - WTP SPACE RANDOM COEFFICIENT MODEL - MURKY WATER DAYS FIRS
> T */
> 
> eststo: mixlogitwtp ce dummy `demos' if block_f == 0 [pw=wgt_final
> ], group(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p m
> d_g md_p) nrep(250) vce(robust) ;

Iteration 0:   log pseudolikelihood =  -1828.034  (not concave)
Iteration 1:   log pseudolikelihood = -1691.5209  (not concave)
Iteration 2:   log pseudolikelihood =  -1679.836  (not concave)
Iteration 3:   log pseudolikelihood = -1678.8979  (not concave)
Iteration 4:   log pseudolikelihood = -1627.6731  (not concave)
Iteration 5:   log pseudolikelihood = -1626.5767  
Iteration 6:   log pseudolikelihood =  -1620.733  
Iteration 7:   log pseudolikelihood = -1620.4563  
Iteration 8:   log pseudolikelihood = -1620.3699  
Iteration 9:   log pseudolikelihood = -1620.3692  
Iteration 10:  log pseudolikelihood = -1620.3692  

Mixed logit model in WTP space                  Number of obs     = 
>      7,496
                                                Wald chi2(12)     = 
>     646.30
Log pseudolikelihood = -1620.3692               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   11.98682   4.571571     2.62   0.009     3.026706  
>   20.94693
        male |   -1.58509   2.302572    -0.69   0.491    -6.098048  
>   2.927868
         age |  -.1321496   .0737773    -1.79   0.073    -.2767505  
>   .0124512
     college |   2.825917   2.268585     1.25   0.213    -1.620428  
>   7.272261
       white |  -3.543359   2.600808    -1.36   0.173     -8.64085  
>   1.554131
        ec_g |   .3472355   .0736919     4.71   0.000      .202802  
>    .491669
        ec_p |  -.3294217   .0768593    -4.29   0.000    -.4800633  
>  -.1787802
        hr_g |   .2212261   .0914891     2.42   0.016     .0419107  
>   .4005416
        hr_p |   -.837515   .2230542    -3.75   0.000    -1.274693  
>  -.4003367
        md_g |   .1151786    .050225     2.29   0.022     .0167393  
>   .2136178
        md_p |   -.293567   .1034775    -2.84   0.005    -.4963791  
>  -.0907549
       mcost |  -1.886546   .1698786   -11.11   0.000    -2.219502  
>   -1.55359
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   .6629562   .0993041     6.68   0.000     .4683238  
>   .8575886
        ec_p |   .1272512   .1380809     0.92   0.357    -.1433824  
>   .3978849
        hr_g |   .6762395    .211125     3.20   0.001      .262442  
>   1.090037
        hr_p |   2.236494    .625222     3.58   0.000     1.011081  
>   3.461906
        md_g |   .5252772   .1639693     3.20   0.001     .2039032  
>   .8466512
        md_p |   .8025447   .1411016     5.69   0.000     .5259906  
>   1.079099
       mcost |  -.2716389   .2812819    -0.97   0.334    -.8229412  
>   .2796634
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est5 stored)

. local ll5 = e(ll) ;

.   local par5 = e(k) ;

. /* HYPOTHESIS TEST - NULL = ORDERING OF ATTRIBUTES DOES NOT MATTER
>  */
> 
> local ll_ur = `ll4' + `ll5' ;

.  local par_ur = `par4' + `par5' ;

. local ll_r  = `ll2' ;

.          local par_r  = `par2' ;

. local test_stat = 2*(`ll_ur' - `ll_r') ;

. local pval  = chi2tail(`par_ur'-`par_r',`test_stat') ;

. local cval  = invchi2(`par_ur'-`par_r',.95) ;

. di "" ;


. di "HYPOTHESIS TEST - NULL = ORDERING OF ATTRIBUTES DOES NOT MATTE
> R" ;
HYPOTHESIS TEST - NULL = ORDERING OF ATTRIBUTES DOES NOT MATTER

. di "TEST STATISTIC:   " `test_stat' ;
TEST STATISTIC:   29.956379

. di "CRITICAL VALUE:   " `cval' ;
CRITICAL VALUE:   30.143527

. di "P-VALUE (5%):     " `pval' ;
P-VALUE (5%):     .05235625

. di "" ;


.      /* #6 - WTP SPACE FIXED COEFFICIENT MODEL - FIRST CHOICE EXPE
> RIMENT ONLY */
> 
> mixlogitwtp ce dummy ec_g ec_p hr_g hr_p md_g md_p `demos' [pw=wgt
> _final] if order == 1, group(coid) id(resp_id) price(mcost) nrep(3
> ) vce(robust) ;

Iteration 0:   log pseudolikelihood = -939.85486  
Iteration 1:   log pseudolikelihood = -939.79603  
Iteration 2:   log pseudolikelihood = -939.79585  
Iteration 3:   log pseudolikelihood = -939.79585  

Mixed logit model in WTP space                  Number of obs     = 
>      3,792
                                                Wald chi2(12)     = 
>     837.47
Log pseudolikelihood = -939.79585               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   37.69579   7.009131     5.38   0.000     23.95815  
>   51.43343
        ec_g |   .2563473   .1728218     1.48   0.138    -.0823773  
>   .5950718
        ec_p |  -.2014403   .1844597    -1.09   0.275    -.5629747  
>   .1600942
        hr_g |  -.0266884    .138435    -0.19   0.847     -.298016  
>   .2446392
        hr_p |  -1.005458   .4623128    -2.17   0.030    -1.911574  
>  -.0993411
        md_g |   .0214812    .115798     0.19   0.853    -.2054787  
>   .2484411
        md_p |  -.0528614   .1961413    -0.27   0.788    -.4372912  
>   .3315685
        male |  -1.000764   2.977144    -0.34   0.737    -6.835858  
>   4.834331
         age |  -.3908643   .0992665    -3.94   0.000    -.5854231  
>  -.1963054
     college |   9.069115   3.417491     2.65   0.008     2.370956  
>   15.76727
       white |   -8.75713   3.946964    -2.22   0.027    -16.49304  
>  -1.021223
       mcost |  -3.143497    .157224   -19.99   0.000     -3.45165  
>  -2.835343
-------------+------------------------------------------------------
> ----------
SD           |
       mcost |   .1590641   .2294721     0.69   0.488     -.290693  
>   .6088211
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive

. matrix define b0 = e(b) ;

. constraint 1 [SD]_b[mcost] = 0 ;

. matrix list b0 ;

b0[1,13]
          Mean:       Mean:       Mean:       Mean:       Mean:
         dummy        ec_g        ec_p        hr_g        hr_p
y1    37.69579   .25634726  -.20144027  -.02668837  -1.0054575

          Mean:       Mean:       Mean:       Mean:       Mean:
          md_g        md_p        male         age     college
y1   .02148118  -.05286136  -1.0007636  -.39086427    9.069115

          Mean:       Mean:         SD:
         white       mcost       mcost
y1  -8.7571296  -3.1434966   .15906406

. matrix b0[1,13] = 0 ;

. eststo: mixlogitwtp ce dummy ec_g ec_p hr_g hr_p md_g md_p `demos'
>  [pw=wgt_final] if order == 1, group(coid) id(resp_id) price(mcost
> ) nrep(3) vce(robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -940.19269  
Iteration 1:   log pseudolikelihood =  -940.1588  
Iteration 2:   log pseudolikelihood = -940.15873  

Mixed logit model in WTP space                  Number of obs     = 
>      3,792
                                                Wald chi2(12)     = 
>     876.01
Log pseudolikelihood = -940.15873               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   38.00801   6.905186     5.50   0.000     24.47409  
>   51.54192
        ec_g |   .2392266   .1646138     1.45   0.146    -.0834106  
>   .5618638
        ec_p |  -.1873303   .1748903    -1.07   0.284    -.5301089  
>   .1554484
        hr_g |  -.0371015   .1338117    -0.28   0.782    -.2993676  
>   .2251645
        hr_p |  -.9636436   .4400835    -2.19   0.029    -1.826191  
>  -.1010958
        md_g |    .024314   .1139269     0.21   0.831    -.1989786  
>   .2476065
        md_p |  -.0499242   .1931757    -0.26   0.796    -.4285417  
>   .3286933
        male |  -.9201195   2.926123    -0.31   0.753    -6.655216  
>   4.814977
         age |  -.3905914    .098545    -3.96   0.000     -.583736  
>  -.1974468
     college |   8.807266   3.240259     2.72   0.007     2.456474  
>   15.15806
       white |  -8.251366   3.604504    -2.29   0.022    -15.31606  
>  -1.186668
       mcost |   -3.12397   .1433597   -21.79   0.000     -3.40495  
>   -2.84299
-------------+------------------------------------------------------
> ----------
SD           |
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est6 stored)

. local ll6 = e(ll) ;

.   local par6 = e(k) ;

. /* #7, #8, #9 - WTP SPACE RANDOM COEFFICIENT MODEL - COUNTY-SPECIF
> IC MODELS */
> 
> mixlogitwtp ce dummy `demos' [pw=wgt_final], group(coid) id(resp_i
> d) price(mcost) rand(ec_g ec_p hr_g hr_p md_g md_p) nrep(50) vce(r
> obust) ;

Iteration 0:   log pseudolikelihood = -3692.9865  (not concave)
Iteration 1:   log pseudolikelihood = -3428.8006  (not concave)
Iteration 2:   log pseudolikelihood = -3358.2181  
Iteration 3:   log pseudolikelihood = -3286.9047  
Iteration 4:   log pseudolikelihood = -3275.9457  
Iteration 5:   log pseudolikelihood = -3275.9062  
Iteration 6:   log pseudolikelihood = -3266.0963  
Iteration 7:   log pseudolikelihood =  -3265.594  (backed up)
Iteration 8:   log pseudolikelihood = -3260.1585  
Iteration 9:   log pseudolikelihood = -3259.6339  
Iteration 10:  log pseudolikelihood = -3259.6277  
Iteration 11:  log pseudolikelihood = -3259.6277  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    1286.50
Log pseudolikelihood = -3259.6277               Prob > chi2       = 
>     0.0000

--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   13.50473   3.756478     3.60   0.000     6.142166  
>   20.86729
        male |  -.7192157   1.657702    -0.43   0.664    -3.968252  
>    2.52982
         age |   -.143351   .0536342    -2.67   0.008    -.2484721  
>  -.0382299
     college |   4.992412   1.900714     2.63   0.009     1.267081  
>   8.717743
       white |  -5.267764   1.930723    -2.73   0.006    -9.051913  
>  -1.483616
        ec_g |   .3141447   .0584945     5.37   0.000     .1994977  
>   .4287918
        ec_p |  -.2772772   .0620811    -4.47   0.000    -.3989538  
>  -.1556005
        hr_g |   .2998475   .0568859     5.27   0.000     .1883532  
>   .4113418
        hr_p |  -.6103604   .1682107    -3.63   0.000    -.9400474  
>  -.2806735
        md_g |   .1351369   .0379608     3.56   0.000      .060735  
>   .2095388
        md_p |  -.1746648   .0747462    -2.34   0.019    -.3211647  
>  -.0281649
       mcost |  -2.005516   .0945948   -21.20   0.000    -2.190919  
>  -1.820114
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |    .852167   .1141431     7.47   0.000     .6284506  
>   1.075883
        ec_p |   .5330561   .2180058     2.45   0.014     .1057726  
>   .9603396
        hr_g |  -.4186542   .1215668    -3.44   0.001    -.6569207  
>  -.1803876
        hr_p |   2.028518   .3939481     5.15   0.000     1.256394  
>   2.800642
        md_g |   .5623094   .1076076     5.23   0.000     .3514024  
>   .7732165
        md_p |   1.009218   .1421348     7.10   0.000     .7306388  
>   1.287797
       mcost |  -.1611785   .1777031    -0.91   0.364    -.5094702  
>   .1871132
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive

. matrix define b0 = e(b) ;

. constraint 1 [SD]_b[mcost] = 0 ;

. matrix list b0 ;

b0[1,19]
          Mean:       Mean:       Mean:       Mean:       Mean:
         dummy        male         age     college       white
y1   13.504727  -.71921569  -.14335099   4.9924123  -5.2677642

          Mean:       Mean:       Mean:       Mean:       Mean:
          ec_g        ec_p        hr_g        hr_p        md_g
y1   .31414472  -.27727717   .29984749  -.61036042   .13513692

          Mean:       Mean:         SD:         SD:         SD:
          md_p       mcost        ec_g        ec_p        hr_g
y1  -.17466479  -2.0055162   .85216699   .53305612  -.41865418

            SD:         SD:         SD:         SD:
          hr_p        md_g        md_p       mcost
y1   2.0285185   .56230943    1.009218   -.1611785

. matrix b0[1,19] = 0 ;

. mixlogitwtp ce dummy `demos' [pw=wgt_final]                 , grou
> p(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p md_g md_
> p) nrep(250) vce(robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -3269.2996  
Iteration 1:   log pseudolikelihood = -3268.6816  
Iteration 2:   log pseudolikelihood = -3266.5766  
Iteration 3:   log pseudolikelihood = -3266.4614  
Iteration 4:   log pseudolikelihood = -3266.4609  
Iteration 5:   log pseudolikelihood = -3266.4609  

Mixed logit model in WTP space                  Number of obs     = 
>     15,168
                                                Wald chi2(12)     = 
>    1536.74
Log pseudolikelihood = -3266.4609               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   15.43363    3.34778     4.61   0.000     8.872106  
>   21.99516
        male |  -1.413475   1.691545    -0.84   0.403    -4.728842  
>   1.901892
         age |  -.1518262   .0497358    -3.05   0.002    -.2493067  
>  -.0543458
     college |   4.113203   1.753979     2.35   0.019     .6754665  
>   7.550939
       white |  -5.942157   1.955212    -3.04   0.002    -9.774302  
>  -2.110012
        ec_g |    .321788   .0597997     5.38   0.000     .2045828  
>   .4389932
        ec_p |  -.2869955   .0617029    -4.65   0.000     -.407931  
>  -.1660601
        hr_g |   .2534983   .0599521     4.23   0.000     .1359943  
>   .3710023
        hr_p |  -.6804636   .1618155    -4.21   0.000    -.9976161  
>  -.3633111
        md_g |   .1471188   .0371599     3.96   0.000     .0742867  
>    .219951
        md_p |  -.1723757   .0738686    -2.33   0.020    -.3171554  
>  -.0275959
       mcost |   -2.01751   .0725419   -27.81   0.000     -2.15969  
>  -1.875331
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   .8558582   .1259692     6.79   0.000     .6089631  
>   1.102753
        ec_p |   .4805592   .2156395     2.23   0.026     .0579137  
>   .9032048
        hr_g |  -.5252376   .1350493    -3.89   0.000    -.7899294  
>  -.2605458
        hr_p |   1.689818   .3763044     4.49   0.000     .9522753  
>   2.427361
        md_g |   .4760443   .1261125     3.77   0.000     .2288683  
>   .7232203
        md_p |      1.045   .1310908     7.97   0.000     .7880664  
>   1.301933
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive

. local llb = e(ll) ;

.   local parb = e(k) ;

. eststo: mixlogitwtp ce dummy `demos' [pw=wgt_final] if county == "
> W", group(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p 
> md_g md_p) nrep(250) vce(robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -1370.0171  
Iteration 1:   log pseudolikelihood = -1363.0945  
Iteration 2:   log pseudolikelihood = -1361.5946  
Iteration 3:   log pseudolikelihood = -1361.5828  
Iteration 4:   log pseudolikelihood = -1361.5828  

Mixed logit model in WTP space                  Number of obs     = 
>      6,280
                                                Wald chi2(12)     = 
>     496.90
Log pseudolikelihood = -1361.5828               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   11.55546   5.020992     2.30   0.021       1.7145  
>   21.39643
        male |  -3.274586   2.205473    -1.48   0.138    -7.597235  
>   1.048062
         age |  -.0578622    .072008    -0.80   0.422    -.1989952  
>   .0832708
     college |   3.573361   2.394444     1.49   0.136    -1.119664  
>   8.266385
       white |  -3.538953   3.027241    -1.17   0.242    -9.472236  
>   2.394329
        ec_g |   .3082428   .0826388     3.73   0.000     .1462738  
>   .4702118
        ec_p |  -.3157663   .0896287    -3.52   0.000    -.4914353  
>  -.1400972
        hr_g |   .3499748   .0818183     4.28   0.000      .189614  
>   .5103357
        hr_p |  -.4990117   .2208193    -2.26   0.024    -.9318095  
>  -.0662138
        md_g |   .0759043   .0524832     1.45   0.148    -.0269609  
>   .1787695
        md_p |  -.2526526   .1167546    -2.16   0.030    -.4814875  
>  -.0238178
       mcost |   -1.84604   .1412065   -13.07   0.000      -2.1228  
>   -1.56928
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   .9065275   .1337247     6.78   0.000     .6444319  
>   1.168623
        ec_p |   .4308816   .1315321     3.28   0.001     .1730835  
>   .6886797
        hr_g |  -.4513792   .1331407    -3.39   0.001    -.7123302  
>  -.1904283
        hr_p |   2.088172   .5536995     3.77   0.000      1.00294  
>   3.173403
        md_g |   .4264345   .1566274     2.72   0.006     .1194504  
>   .7334186
        md_p |   1.054922   .2890882     3.65   0.000     .4883196  
>   1.621524
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est7 stored)

. local ll7 = e(ll) ;

.   local par7 = e(k) ;

. eststo: mixlogitwtp ce dummy `demos' [pw=wgt_final] if county == "
> G", group(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p 
> md_g md_p) nrep(250) vce(robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -616.81964  
Iteration 1:   log pseudolikelihood = -612.33949  
Iteration 2:   log pseudolikelihood = -609.27307  (not concave)
Iteration 3:   log pseudolikelihood = -608.95992  
Iteration 4:   log pseudolikelihood = -608.80248  
Iteration 5:   log pseudolikelihood = -608.80145  
Iteration 6:   log pseudolikelihood = -608.80145  

Mixed logit model in WTP space                  Number of obs     = 
>      4,720
                                                Wald chi2(12)     = 
>     578.50
Log pseudolikelihood = -608.80145               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   18.33648   5.775468     3.17   0.001     7.016774  
>   29.65619
        male |  -2.150907   2.814129    -0.76   0.445    -7.666498  
>   3.364685
         age |  -.2021053   .0811233    -2.49   0.013     -.361104  
>  -.0431065
     college |   1.223412    2.67747     0.46   0.648    -4.024334  
>   6.471157
       white |  -4.450112   3.243202    -1.37   0.170    -10.80667  
>   1.906447
        ec_g |   .2888091   .0942642     3.06   0.002     .1040546  
>   .4735636
        ec_p |  -.2276054   .0923543    -2.46   0.014    -.4086166  
>  -.0465942
        hr_g |   .2092729   .1104459     1.89   0.058    -.0071971  
>   .4257429
        hr_p |  -.5586183   .2282215    -2.45   0.014    -1.005924  
>  -.1113124
        md_g |   .1033144   .0649462     1.59   0.112    -.0239778  
>   .2306067
        md_p |  -.5009067   .1446837    -3.46   0.001    -.7844816  
>  -.2173319
       mcost |  -2.035582   .1301394   -15.64   0.000    -2.290651  
>  -1.780514
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |    .585545   .2225454     2.63   0.009     .1493641  
>   1.021726
        ec_p |  -.0441274   .2002418    -0.22   0.826    -.4365942  
>   .3483394
        hr_g |  -.7251031   .2458873    -2.95   0.003    -1.207033  
>  -.2431729
        hr_p |   1.879099   .6972976     2.69   0.007     .5124211  
>   3.245777
        md_g |   .3779299   .1739287     2.17   0.030     .0370359  
>    .718824
        md_p |   1.086365   .2221432     4.89   0.000     .6509724  
>   1.521758
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est8 stored)

. local ll8 = e(ll) ;

.   local par8 = e(k) ;

. eststo: mixlogitwtp ce dummy `demos' [pw=wgt_final] if county == "
> M", group(coid) id(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p 
> md_g md_p) nrep(250) vce(robust) constraints(1) from(b0) ;

Iteration 0:   log pseudolikelihood = -1278.4722  
Iteration 1:   log pseudolikelihood = -1265.5882  
Iteration 2:   log pseudolikelihood =  -1262.916  
Iteration 3:   log pseudolikelihood = -1262.8719  
Iteration 4:   log pseudolikelihood = -1262.8718  

Mixed logit model in WTP space                  Number of obs     = 
>      4,168
                                                Wald chi2(12)     = 
>     678.92
Log pseudolikelihood = -1262.8718               Prob > chi2       = 
>     0.0000

 ( 1)  [SD]mcost = 0
--------------------------------------------------------------------
> ----------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
>  Interval]
-------------+------------------------------------------------------
> ----------
Mean         |
       dummy |   13.48452    6.21876     2.17   0.030     1.295975  
>   25.67307
        male |   2.504419   3.042373     0.82   0.410    -3.458523  
>   8.467361
         age |  -.2050789   .0894086    -2.29   0.022    -.3803166  
>  -.0298411
     college |   7.199573   3.480634     2.07   0.039     .3776568  
>   14.02149
       white |  -7.640246   3.950606    -1.93   0.053    -15.38329  
>   .1027996
        ec_g |   .3484245   .1216916     2.86   0.004     .1099134  
>   .5869357
        ec_p |   -.237542   .1119291    -2.12   0.034    -.4569191  
>  -.0181649
        hr_g |   .2515502   .1205164     2.09   0.037     .0153424  
>   .4877581
        hr_p |  -.9452335    .324374    -2.91   0.004    -1.580995  
>  -.3094722
        md_g |   .2865357   .0851612     3.36   0.001     .1196227  
>   .4534486
        md_p |   .1167976   .1389382     0.84   0.401    -.1555163  
>   .3891114
       mcost |  -2.115651   .1190652   -17.77   0.000    -2.349015  
>  -1.882288
-------------+------------------------------------------------------
> ----------
SD           |
        ec_g |   1.011927   .1862751     5.43   0.000     .6468347  
>    1.37702
        ec_p |   .4443315   .3463257     1.28   0.199    -.2344543  
>   1.123117
        hr_g |  -.3534234   .2481117    -1.42   0.154    -.8397135  
>   .1328667
        hr_p |   1.215197   .5820369     2.09   0.037     .0744255  
>   2.355968
        md_g |   .8175701   .3043679     2.69   0.007       .22102  
>    1.41412
        md_p |   .8226355   .3041111     2.71   0.007     .2265887  
>   1.418682
       mcost |          0  (omitted)
--------------------------------------------------------------------
> ----------
The sign of the estimated standard deviations is irrelevant: interpr
> et them as
being positive
(est9 stored)

. local ll9 = e(ll) ;

.   local par9 = e(k) ;

. /* HYPOTHESIS TEST - NULL = NO PREFERENCE HETEROGENEITY ACROSS COU
> NTIES */
> 
> local ll_ur = `ll7' + `ll8' + `ll9' ;

.  local par_ur = `par7' + `par8' + `par9' - 3 ;

. local ll_r  = `llb' ;

.          local par_r  = `parb' - 1 ;

. local test_stat = 2*(`ll_ur' - `ll_r') ;

. local pval  = chi2tail(`par_ur'-`par_r',`test_stat') ;

. local cval  = invchi2(`par_ur'-`par_r',.95) ;

. di "" ;


. di "HYPOTHESIS TEST - NULL = NO PREFERENCE HETEROGENEITY ACROSS CO
> UNTIES" ;
HYPOTHESIS TEST - NULL = NO PREFERENCE HETEROGENEITY ACROSS COUNTIES

. di "TEST STATISTIC:   " `test_stat' ;
TEST STATISTIC:   66.409606

. di "CRITICAL VALUE:   " `cval' ;
CRITICAL VALUE:   50.99846

. di "P-VALUE (5%):     " `pval' ;
P-VALUE (5%):     .00150543

. di "" ;


.      /* #10 - NOW MODEL WITH ALL RESPONDENTS */
> 
> use $ce_data, replace ;

. /* DROPPING:                                                      
>                       */
> /*    1) SPEEDING PEOPLE (<=8 MINUTES TO COMPLETE SURVEY)         
>                       */
> /*    2) PEOPLE WHO TAKE MORE THAN A WEEK TO COMPLETE THE SURVEY  
>                       */
> /*    3) PEOPLE WHO FAIL THE TRAP QUESTION                        
>                       */
> /*    4) PEOPLE WHO DISAGREE OR STRONGLY DISAGREE WITH THE CONSEQU
> ENTIALITY STATEMENTS  */
> 
> *keep if survey_duration > 8 & survey_duration <= (7*24*60) & inli
> st(debrief_3,"D") & !inlist(debrief_1,"D","SD") & !inlist(debrief_
> 2,"D","SD") ;
. rename gender male ;

. local demos = "male age college white" ;

. sort resp_id seq ;

. preserve ;

.    by resp_id : keep if _n == 1 ;
(7,327 observations deleted)

.    log on ;
(log already on)

.       sum `demos' ;

    Variable |        Obs        Mean    Std. Dev.       Min        
> Max
-------------+------------------------------------------------------
> ---
        male |      2,511    .4752197    .4838155          0        
>   1
         age |      2,511    51.14356    16.14859         18        
>  90
     college |      2,511    .7237811    .4338343          0        
>   1
       white |      2,511    .7395972    .4254164          0        
>   1

.    log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Fin
> al_Data_Public\log\5-sensitivities.log
  log type:  text
 paused on:  21 Jan 2023, 21:45:14
--------------------------------------------------------------------
--------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Fin
> al_Data_Public\log\5-sensitivities.log
  log type:  text
resumed on:  21 Jan 2023, 21:45:14

. foreach y of varlist `demos' { ;
  2.    replace `y' = 0 if dummy == 0 ;
  3. } ;
(4,894 real changes made)
(9,838 real changes made)
(7,230 real changes made)
(7,418 real changes made)

. log on ;
(log already on)

. /* #10 - WTP SPACE RANDOM COEFFICIENT MODEL - ALL RESPONDENTS */
> 
> eststo: mixlogitwtp ce dummy `demos' [pw=wgt_final], group(coid) i
> d(resp_id) price(mcost) rand(ec_g ec_p hr_g hr_p md_g md_p) nrep(2
> 50) vce(robust) ;

Iteration 0:   log pseudolikelihood = -5159.9561  (not concave)
Iteration 1:   log pseudolikelihood = -4637.8414  (not concave)
Iteration 2:   log pseudolikelihood = -4527.8071  
Iteration 3:   log pseudolikelihood = -4480.9389  (not concave)
Iteration 4:   log pseudolikelihood = -4477.5954  
Iteration 5:   log pseudolikelihood = -4474.1161  
Iteration 6:   log pseudolikelihood = -4474.0198  (n
> ot concave)
Iteration 7:   log pseudolikelihood = -4473.9874  (not concave)
Iteration 8:   log pseudolikelihood = -4473.1803  (not concave)
Iteration 9:   log pseudolikelihood = -4472.2979  
Iteration 10:  log pseudolikelihood = -4470.7847  
Iteration 11:  log pseudolikelihood =  -4470.762  
Iteration 12:  log pseudolikelihood =  -4470.762  

Mixed logit model in WTP space                  Number of obs     =     19,676
                                                Wald chi2(12)     =    1126.20
Log pseudolikelihood =  -4470.762               Prob > chi2       =     0.0000

------------------------------------------------------------------------------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Mean         |
       dummy |   19.15369   2.931296     6.53   0.000     13.40846    24.89893
        male |   -3.51306   1.553498    -2.26   0.024    -6.557859   -.4682601
         age |  -.1661274      .0435    -3.82   0.000    -.2513859   -.0808689
     college |   2.535878    1.49984     1.69   0.091    -.4037551    5.475511
       white |  -6.003408   1.647562    -3.64   0.000     -9.23257   -2.774246
        ec_g |   .2515157   .0509185     4.94   0.000     .1517173    .3513141
        ec_p |  -.2644068   .0519118    -5.09   0.000    -.3661521   -.1626616
        hr_g |   .2165678   .0550966     3.93   0.000     .1085804    .3245552
        hr_p |  -.4923037   .1359301    -3.62   0.000    -.7587218   -.2258855
        md_g |   .1283642   .0312468     4.11   0.000     .0671217    .1896068
        md_p |  -.1349586   .0667045    -2.02   0.043    -.2656969   -.0042202
       mcost |  -1.866928   .1191614   -15.67   0.000     -2.10048   -1.633376
-------------+----------------------------------------------------------------
SD           |
        ec_g |   .8150723   .0724016    11.26   0.000     .6731677    .9569768
        ec_p |   .3463597   .1117755     3.10   0.002     .1272837    .5654356
        hr_g |   .5901878   .0871586     6.77   0.000       .41936    .7610156
        hr_p |   2.131738   .2472244     8.62   0.000     1.647187    2.616289
        md_g |   .4425091    .127095     3.48   0.000     .1934074    .6916107
        md_p |   1.123512   .1156702     9.71   0.000     .8968027    1.350222
       mcost |    .423159   .1633366     2.59   0.010     .1030252    .7432928
------------------------------------------------------------------------------
The sign of the estimated standard deviations is irrelevant: interpret them as
being positive
(est10 stored)

. local ll10 = e(ll) ;

.   local par10 = e(k) ;

. esttab using log\5-sensitivities.rtf, replace title(Alternative Specifications) mtitles("WTP-SP
> ACE - FIXED" "WTP-SPACE - RANDOM" "NO WEIGHTS" "ECO COND FIRST" "MURKY WATER FIRST" "FIRST CE O
> NLY" "WAKE CO ONLY" "MECKLENBERG CO ONLY" "GUILFORD CO ONLY" "ALL RESPONDENTS") scalars(ll) sta
> r(* .10 ** .05 *** .01) ;
(output written to log\5-sensitivities.rtf)

. log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\5-sensitiv
> ities.log
  log type:  text
 paused on:  21 Jan 2023, 23:29:18
-------------------------------------------------------------------------------------------------
